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ABSTRACT 

We present analyses of the ASCA X-ray spectra of two Seyfert galaxies, Tololo 0109- 
383 and ESO 138-G1. In both cases, spectral fitting reveals two statistically acceptable 
continuum models: Compton reflection and partial covering. Both spectra have strong 
iron Ka lines, with equivalent widths greater than 1.5 keV. These large equivalent 
widths are suggestive of heavier obscuration than that directly indicated by the partial- 
covering models (« 2 x 10 23 cm -2 ), with the actual column densities being 'Compton- 
thick' (i.e. TVh <^ 1.5 x 10 24 cm -2 ). We use the hard X-ray/[0 III] flux correlation for 
Seyferts and data from the literature to provide additional support for this hypothesis. 
Since Tololo 0109-383 is known to have optical type 1 characteristics such as broad 
Balmer line components and Fe n emission, this result marks it as a notable object. 

Key words: galaxies: individual: Tololo 0109-383 - galaxies: individual: ESO 138- 
Gl - galaxies: nuclei - galaxies: Seyfert - X-rays: galaxies. 



1 INTRODUCTION 

The issue of absorption in Seyfert nuclei is one that has 
motivated a great deal of research. For example, studying 
the absorbing gas in Seyfert 2s offers the most straightfor- 
ward means for learning about the putative molecular torus 
surrounding the central engine, which is a cornerstone of 
the unified model for Seyferts (e.g. Antonucci 1993). These 
studies are also pertinent to research aiming to discover the 
source of the cosmic X-ray background radiation. In both 
of these cases, a key issue is to understand the distribution 
of absorbing column densities of the Seyfert population. A 
recent study by Risaliti, Maiolino & Salvati (1999) reported 
that roughly half of Seyfert 2s have 'Compton-thick' intrin- 
sic absorption columns (-/Vh ii c^ 1 = 1.5 x 10 24 cm -2 ), a 
fraction that is higher than indicated by earlier studies that 
were biased toward bright X-ray sources. Since the number 
of known and well-studied Compton-thick Seyferts is not 
large, the study of new nearby examples is significant. In 
this paper, we present a study of two Seyferts, Tololo 0109- 
383 (NGC 424) and ESO 138-G1, in which we find evi- 
dence for the presence of Compton-thick nuclear absorp- 
tion. ESO 138-G1 is a Seyfert 2 (e.g. Alloin et al. 1992). 
Tololo 0109-383 was originally classified as a Seyfert 2 as 
well (Smith 1975), but later studies (Boisson & Durret 1986; 
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Table 1. Basic source properties. 



Object 


V 


2 


Galactic jVjj 


Tololo 0109-383 


13.9 


0.012 


1.8 X 10 20 cm" 2 a 


ESO 138-G1 


14.3 


0.0091 


1.6 X 10 21 cm" 2 b 



a Stark et al. (1992) b Heiles & Cleary (1979) 



Durret & Bergeron 1988; Murayama, Taniguchi & Iwasawa 
1998) cast this into doubt. Based on the presence of a broad 
Balmer line component and Fe n emission lines in its optical 
spectrum, Murayama et al. (1998) argued for the type 1 na- 
ture of Tololo 0109-383, settling on a marginal classification 
between type 1 and type 2. One of our goals in this study 
is to use its X-ray properties to gain insight into its actual 
nature. 

Both Tololo 0109-383 and ESO 138-G1 are relatively 
nearby (see Table 1) and display the so-called coronal lines, 
high-ionization forbidden optical emission lines from species 
such as [Fe vn], [Fe x], and [Fe xiv] (Alloin et al. 1992; Mu- 
rayama et al. 1998). Tololo 0109-383 is interesting in that 
it is one of only a few known galaxies in which the coro- 
nal line region has been observed to be spatially extended. 
In this study, however, we concern ourselves strictly with 
these galaxies' X-ray emission and their optical characteris- 
tics (such as [O ill] A5007 emission) as they relate to their 
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Energy (keV) 

Figure 1. Plots of Model A data-to-model ratio for (a) 
Tololo 0109-383 and (b) ESO 138-G1. Note the strong residu- 
als due to iron Ka emission from 6-7 keV. Dots are SIS data 
points, and plain crosses arc GIS data points. 



X-ray properties. We adopt Ho = 70 km s Mpc and 
go = §. 



2 X-RAY OBSERVATIONS AND ANALYSIS 

2.1 Observations and data reduction 

The ASCA observations of Tololo 0109-383 and ESO 138- 
Gl were performed on 2-3 July 1997 and 6-7 Septem- 
ber 1997, respectively. For our analysis, we used Revision 2 
processed data from Goddard Space Flight Center, prepared 
using standard screening criteria (Pier 1997). After screen- 
ing, the respective SIS/GIS exposure times were 34 ks/34 ks 
for Tololo 0109-383 and 27 ks/31 ks for ESO 138-G1. 
For both objects, we utilized ximage (Giommi, Angelini & 
White 1997) to locate the sources in the SIS and GIS im- 
ages. We made use of xselect (Ingham & Arnaud 1998) to 
reduce the data. 



2.2 Spectral analysis 

In order to permit \ 2 fitting, we adopted a minimum 
spectral group size of 15 events per data point. For the 
SIS/GIS detectors, we have included the energy ranges 0.6- 
10 keV/0.9-10 keV. We performed simultaneous fitting of 
the spectra from all four detectors for each object, and we 
constrained all parameters to be the same for each group, 
aside from the absolute model normalizations. We report all 
equivalent widths and fluxes from the SIS0 detector, and 
all errors at the 90% confidence level unless otherwise indi- 
cated. We have used xspec (Arnaud 1996) for all spectral 
fitting. See Table 2 for a summary of the \ 2 fitting. 

We began by fitting the spectra of both objects with 
a simple power law absorbed by the Galactic column den- 
sity (Model A). We have allowed for uncertainties of up to 
20% in the Galactic column densities; such variations do 
not materially alter our results. Model A fails to provide a 
statistically acceptable fit for either spectrum. As shown in 
Figure 1, there are strong line- like residuals above 6 keV 
in both spectra. In addition, the photon indices are sig- 
nificantly flatter than expected for intrinsic Seyfert spectra 
(0.99tg:i4 for Tololo 0109-383 and 0.551^3 for ESO 138- 
Gl). In Model B, we constrained the photon index of the 
power law to lie in the canonical range for Seyferts; specif- 
ically, we required it to be above 1.6 (e.g. Brandt, Mathur 
& Elvis 1997). We also added a Gaussian emission line to 
represent iron Ka emission. We constrained the line width 
parameter a to lie below 0.4 keV in order to remain consis- 
tent with known Seyfert characteristics (e.g. Nandra et al. 
1997). The fit improves somewhat for Tololo 0109-383 and 
is comparable to Model A for ESO 138-G1 (see Table 2); 
neither is statistically acceptable. Our next step was to add 
an intrinsic absorption column (Model C), but the best-fit 
value of this column density is zero for Tololo 0109-383, and 
therefore the fit does not improve. The improvement in the 
fit of the spectrum for ESO 138-G1 is slight. 

We were able to achieve statistically acceptable fits of 
both spectra using models consisting of absorption by the 
Galactic column densities of Compton-refiected continua 
('pexrav') plus iron lines (Model D). This type of model 
is often found to represent well the spectra of absorbed 
Seyferts. Because of the poor signal-to-noise of our data, 
we constrained the element abundances to be solar, and we 
fixed the power-law cutoff energy to be far outside the ASCA 
band at 1000 keV. In Table 3 we list the relevant fit parame- 
ters. The high values of the reflection scaling factors indicate 
that reflected components dominate both spectra, and the 
photon indices, although poorly constrained, are consistent 
with reasonable values. We note that the w 6.4 keV lines 
are consistent with neutral iron Ka emission and that the 
statistical degradation of the fit is not large if we constrain 
the line in the spectrum of Tololo 0109-383 to be narrow 
(a = 0.05 keV). 

Statistically acceptable fits can also be achieved using 
models consisting of power laws plus iron lines absorbed by 
the Galactic column densities and partial-covering (intrin- 
sic) absorption columns (Model E). Like reflection models, 
partial covering is also often found to be important in the 
spectra of Seyferts with obscuration. The photon indices and 
line equivalent widths obtained from this model are consis- 
tent with those obtained from Model D, and again the line 
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Table 2. Model statistics. 



Object Name 


Parameter 


Model A 


Model B 


Model C 


Model D 


Model E 


Model F 


Tololo 0109-383 


X 2 /dof 


288 / 171 


208 / 171 


208 / 171 


164 / 171 


158 / 171 


190 / 171 




1 - P(x 2 | vT 


5.0 x 10~ 13 


9.5 x 10~ 3 


8.2 x 10~ 3 


0.44 


0.54 


5.0 x 10~ 2 


ESO 138-G1 


X 2 /dof 


434 / 289 


434 / 289 


412 / 289 


288 / 289 


281 / 289 


376 / 289 




1 - P( x 2 | v) 


8.3 x lO" 30 


1.2 x 10~ 8 


4.4 x 10~ 7 


0.34 


0.46 


6.8 x 10~ 5 



a P(x 2 I v) is the probability that the observed chi-square for a correct model should be less than a value x 2 (Press 
ct al. 1989). 



Table 3. Best-fit Model D parameters with 90% (Ax 2 
errors, X-ray fluxes, and [O III] fluxes. 



2.71) 



Parameter 21 

Refl. scaling factor 
Photon index (r) 
Line energy (kcV) 
Line width (<r) (keV) 
Equivalent width (keV) 
F h 

^2-10 

F c 
[O in] 



Tololo 0109-383 ESO 138-G1 



2.3 
6.43 
0.20 

1. 



+ 120 
-14 
+0.2 
0.2 
+0.09 
0.10 
+0.28 
0.12 
+0.4 



70: 



-0.5 

2.0 x 10- 

1.2 x 10" 

2.3 x 10- 



130 
32 
2 1+ ' 4 

6 36+ 05 
07+ - 09 
1 7 +o.4 

1 -'-0.4 

1.8 x 10~ 13 
1.8 x 10~ 12 
9.3 x 10~ 13 



a All fluxes are reported in erg cm -2 s _1 . b SIS0 observed 0.6- 
2 keV and 2—10 keV fluxes. c [0 III] A5007 emission, corrected for 
Galactic extinction. From Murayama et al. (1998) and Schmitt 
& Storchi-Bcrgmann (1995), respectively. The value reported for 
ESO 138-G1 is the total flux within 5" of the nucleus. 



energies and widths are consistent with neutral, narrow iron 
Ka emission. We shall later discriminate between Models D 
and E based on physical arguments. 

As a means of checking the general robustness of these 
results, we also considered alternate models consisting of 
Raymond-Smith plasma components added to Model C 
(Model F). Physically the plasma component represents 
thermal gas emission such as is often associated with star- 
burst activity (e.g. Ptak et al. 1999). For each spectrum, the 
fit is a significant improvement over Model B; however, the 
fits are still not statistically acceptable. 



3 DISCUSSION AND CONCLUSIONS 

3.1 The iron Ka lines and variability 

Iron Ka lines are produced through the reprocessing of 
primary X-rays and become strongest in equivalent width 
when the primary X-ray continuum is suppressed in the 
neighborhood of 6.4-6.97 keV. The large iron Kq equivalent 
widths (greater than 1.5 keV) and flat apparent spectral con- 
tinua (see Section 2.2) that we observe in Tololo 0109-383 
and ESO 138-G1 indicate that reprocessing is important in 
these sources (e.g. Matt, Brandt & Fabian 1996), a deduc- 
tion which seems to be confirmed by the good fits of the 
spectra that we achieved using Compton-reflection models 
(Model D). The large reflection scaling factors in those fits 
indicate that the reflected components dominate any direct 
components of the nuclear emission to a great extent. We 
achieved equally good fits, however, using partial-covering 
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Figure 2. ASCA SIS (dots) and GIS (plain crosses) spectra of 
(a) Tololo 0109-383 and (b) ESO 138-G1. Model D has been fit 
to the data. The ordinates for the lower panels (labeled x) show 
the fit residuals in terms of sigmas with error bars of size one. 



models. We must therefore discriminate between the two 
models based on physical arguments. 

Matt et al. (1996) describe two reprocessing mecha- 
nisms that can be significant in Seyferts, both of which can 
produce very strong iron Ka emission. In the first case, the 
iron emission occurs at the inner surface of the circumnu- 
clear absorbing material (the putative torus). The viewing 
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Figure 3. F2-10 versus F^q nI j for the Seyferts in the Mulchacy 
ct al. (1994) sample plus Tololo 0109-383 and ESO 138-G1. We 
have updated the 2-10 kcV flux for NGC 4051 (a Seyfert 1) using 
the average flux reported by Guainazzi ct al. (1996). The object in 
the lower right-hand corner, NGC 1068, is a well-known Compton- 
thick Seyfert 2 (e.g. Matt ct al. 1997). 



angle can be such that the central source is obscured but the 
inner wall of the torus is not. In this case, predominantly 
neutral emission is expected and iron Kq equivalent widths 
can be up to a few keV. Our Model D is entirely consistent 
with this theoretical mechanism. The second scenario places 
the iron emission in optically thin gas in the nuclear region 
(the 'warm mirror'). As long as the direct X-ray emission is 
obscured, the iron Kq equivalent width can be as high as a 
few keV, but in this case the line is expected at higher ener- 
gies (6.7-6.97 keV), representing emission from very highly 
ionized iron. If we interpret the partial-covering column den- 
sities as representing absorbing gas along some lines of sight 
through the scattering regions, our Model E is consistent 
with this scenario, except for the energies of the iron Ka 
lines. These energies are inconsistent at greater than the 
90% confidence level with emission from the highest ioniza- 
tion states of iron. Thus we can conclude that a scattered 
component is not highly significant in the spectra of either 
Tololo 0109-383 or ESO 138-G1. 

Based on these considerations, it is likely that we are 
primarily seeing reprocessing by the inner surface of a cir- 
cumnuclear absorber (torus). In order to get the large equiv- 
alent widths we observe (greater than 1.5 keV), however, 
an absorption column density is required that is signifi- 
cantly larger than the intrinsic column densities indicated 
by Model E. A column density of 2 x 10 23 cm~ 2 (of the 
order we fit in Tololo 0109-383 and ESO 138-G1) absorbs 
only « 25% of the X-ray flux at 6.4 keV. To get such a large 
equivalent width in the 6-7 keV range, the intrinsic absorp- 
tion along a direct line of sight to the X-ray source must be 
essentially Compton-thick (e.g. Matt et al. 1996). 

We have also analyzed ROSAT data for Tololo 0109- 
383 and ESO 138-G1 in order to check for flux variability 



between the ROSAT and ASCA observations of these ob- 
jects. Tololo 0109-383 was observed by the ROSAT PSPC 
on 2-3 July 1992 (Rush & Malkan 1996), and ESO 138-G1 
appears in the ROSAT HRI field of SN1990W, which was 
observed on 11-12 September 1995. Hence the variability 
timescales we are probing are roughly 5 years and 2 years, 
respectively. Significant low-energy variability would place 
a constraint upon the size of the scattering/reflecting re- 
gions. Based on our analysis of the data and the Fo.1-2 for 
Tololo 0109-383 quoted by Rush & Malkan (1996), how- 
ever, we find no evidence for significant variability between 
the ROSAT and ASCA observations of either source. The 
fluxes are consistent to within the cross-calibration uncer- 
tainties of the detectors (e.g. Iwasawa, Fabian & Nandra 
1999). 



3.2 Additional evidence for Compton-thick 
absorption 

We now consider additional evidence supporting the conclu- 
sion that the primary nuclear sources in Tololo 0109-383 and 
ESO 138-G1 are completely obscured in the ASCA band. 
Mulchaey et al. (1994) reported correlations between Fy Q ln ] 
(see Table 3) and F2-10 for Seyferts. They found the follow- 
ing relationships (the reported error ranges are one standard 
deviation) : 



l0g(J[O m]/f2- 
l0g(F[O m\/F2- 



10) = -1.89 ±0.50 
10) = -1.76 ±0.62 



(Seyfert Is) 
(Seyfert 2s) 



A simple calculation using the fluxes from Table 3 yields the 
values log(F [0 mj/i^-io) = -0.72 for Tololo 0109-383 and 
log(F [0 ln] /F2-io) = -0.29 for ESO 138-G1. These num- 
bers fall well outside the la and 2a ranges, respectively, of 
even the Seyfert 2 relationship (see Figure 3). This provides 
further support for the obscuration of most of the direct 
X-rays from these sources below 10 keV, strengthening the 
conclusion that the intrinsic absorption in these galaxies is 
Compton-thick. To emphasize this point, we have used the 
Seyfert 2 relationship to predict F2-10 for the two galaxies 
from F[ m] (e.g. Turner et al. 1997), and we compare these 
to the observed values of Fb-io- Using the Table 3 fluxes, 
we obtain values of 11+f for Tololo 0109-383 and 30+23 for 
ESO 138-G1 for the ratios of predicted to observed fluxes. 
We note that this is a somewhat conservative estimate for 
Tololo 0109-383, to which we might also justifiably apply 
the Seyfert 1 relationship. 



3.3 Tololo 0109-383, a Compton-thick type 1 
Seyfert 

Although the optical properties of Tololo 0109-383 support 
a Seyfert 1 or intermediate nature (see Section 1), its in- 
ferred absorption column density and iron Kq line are more 
like those of Seyfert 2s. It seems clear that in the case of this 
object, the simplified picture of a type 1 or type 2 nature 
falls somewhat short of the truth. Nevertheless, we argue 
that Tololo 0109-383 is a Compton-thick Seyfert; this seems 
to be the only plausible explanation for its observed X-ray 
properties. The combination of its type 1 optical properties 
with this heavy intrinsic absorption marks it as a member 
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of a peculiar class of objects. It may well serve as a nearby 
archetype for the heavily obscured type 1 objects recently 
found in deep X-ray surveys (e.g. Comastri et al. 2000; also 
see Brandt, Laor & Wills 2000). As we have only placed a 
lower limit on the intrinsic absorption column in this ob- 
ject, we point out that the remaining uncertainty in its ac- 
tual value should be resolvable by an observation further 
into the hard X-ray band by a satellite such as BeppoSAX . 
Finally, we comment that the existence of type 1 objects 
with heavy X-ray obscuration, such as Tololo 0109-383 and 
Broad Absorption Line QSOs (e.g. Gallagher et al. 1999), 
should give pause to those who would classify an object as 
type 2 based solely on a hard X-ray spectrum. 
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